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Abstract 



We derive a normal form for a near-integrable, four-dimensional symplectic map with a 
fold or cusp singularity in its frequency mapping. The normal form is obtained for when 
the frequency is near a resonance and the mapping is approximately given by the time-T 
mapping of a two-degree-of freedom Hamiltonian flow. Consequently there is an energy-like 
invariant. The fold Hamiltonian is similar to the well-studied, one-degree-of freedom case 
but is essentially nonintegrable when the direction of the singular curve in action does not 
coincide with curves of the resonance module. We show that many familiar features, such 
as multiple island chains and reconnecting invariant manifolds, are retained even in this 
case. The cusp Hamiltonian has an essential coupling between its two degrees of freedom 
even when the singular set is aligned with the resonance module. Using averaging, we 
approximately reduced this case to one degree of freedom as well. The resulting Hamiltonian 
and its perturbation with small cusp-angle is analyzed in detail. 

1 Introduction 

In a Liouville integrable Hamiltonian flow with d degrees of freedom almost all motion takes 
place on invariant d-dimensional tori. The motion on these tori is conjugate to a linear flow 
with frequencies uj G M°'. The linear flow is {Oo,t) ^ 6 = ojt + 9q where G T'^ is an angle. The 
map that assigns to each torus its frequency lo is called the frequency map. For a Hamiltonian 
flow there are two natural frequency maps. In addition to the one just defined, the frequency- 
ratio map is obtained by fixing the energy and by turning the flow into a Poincare first return 
map. The resulting map has n = d — 1 degrees of freedom, and its frequencies are the n 
frequency ratios loi : uj2 of the original flow. In this way the frequency map of an 

integrable map can always be interpreted as the frequency-ratio map of some flow. Typically 
these are anharmonic — the oscillation frequency changes with the amplitude of oscillation. In 
Hamiltonian dynamics this property is called "twist." When this twist property breaks down, we 
say the frequency map has a singularity. Twist singularities are are not unusual; for example, 
the frequency is not a monotone function of action in any system that has a pair of nested 
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separatrices |M.B991 IMor02j . This "fold" singularity also generically occurs in any system near 
tripling resonances |Moe90l IDMS991 IDM03j . In this paper we will discuss both the fold and 
cusp singularities and some of the dynamical consequences of the breakdown of twist. The cusp 
singularity corresponds to the collision of a pair of fold curves; it is another stable singularity 
when n > 2. 

Consider the integrable symplectic map {9, J) [6 + i}{J),J) with € T" and J € M". 
Because the map is symplectic, the frequency map is Lagrangian, which means that is the 
gradient of a scalar generating function S{J). The twist corresponds to the Jacobian matrix 
r = DO,, or equivalently to the Hessian D^S. There is a singularity in the frequency map 
when det(r) = for some set of action values. Such a set is called the singular set and the 
corresponding tori are said to have vanishing twist. The image of the singular set is the caustic. 
A fold is one such singularity and it is stable for Lagrangian mappings; consequently, a one 
parameter family of such maps will not destroy the fold, but move it around in the frequency 
space. In particular, the caustic will cross rational frequencies. When the map is perturbed by 
a small periodic perturbation interesting dynamics is expected when the fold is near a resonant 
frequency such that the resonance is in the image of the frequency map. This situation is well 
studied for the case n = 1 where the map is area-preserving and leads to a simple one-degree- 
of freedom Hamiltonian model [IIH95, Sim9§j. The canonical example is called the "standard 
nontwist map," for reviews see (AWM03i iWAM041 IWAFM05] . In our study we will derive a 
generalization of this map for the four dimensional case. 

The transformation of the map to a normal form is inspired by singularity theory of La- 
grangian mappings |A(;ZV851 IAKN88lfI7Tl95] . The class of transformat ions is, however, more 
restricted. The canonical variables of an integrable map are the angle-action coordinates so the 
transformations should respect the periodicity of the angles. Moreover, since we are in a near- 
integrable situation there is a natural distinction between momenta and coordinates which we 
would like to maintain as much as possible. These restrictions are quite severe, for example they 
restrict the choice of coordinate transformations to cotangent lifts of linear maps A G SL(n,7j). 
In singularity theory, linear transformations are used to normalize the quadratic terms, while 
nonlinear transformations are used to remove most of the cubic terms. The latter is impossible 
in our setting, since we are restricted to linear transformations, but we can use scaling to bring 
out the singular behavior. Moreover, we cannot diagonalize the quadratic part using orthogonal 
transformations, since these are in general not in 5'L(n,Z). 

This will affect the normal form that we derive when the singular curve is not aligned with 
the resonance: we find it better to use noncanonical variables to decouple the required angle 
and action transformations. The resulting normal form is a Hamiltonian system, but one whose 
Poisson bracket is not in standard form. 

After a brief review of the fold and cusp singularities in the context of integrable symplectic 
maps in ^ we will derive the normal form for a weakly nonintegrable mapping in ^ The 
resulting normal form has a two degree-of-freedom Hamiltonian, shown in ^ The dynamical 
behavior for the fold case and the cusp case are studied in the last two sections, ^and ^ 

2 Fold and Cusp Singularities for Integrable Symplectic Maps 

An integrable, four-dimensional, symplectic map in angle-action coordinates has the form 
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where G x M^. 

The dynamics of this system are completely encoded in the frequency map $7 : — > R^, 
that defines the rotation vector on each invariant torus, I = c. This map is symplectic with the 
standard form w = dli A difi + d/2 A dip2 providing that the Jacobian matrix DO. is symmetric. 

The map is determined by a mixed-variable generating function 



through the implicit equations 99' = dS/dl' and / = dS/dip. Therefore, the gradient of the 
"kinetic energy" T : — > M gives the the frequency map, 



Note that this implies that the Jacobian of Q is symmetric since it is the Hessian of T. 

In this paper we will study perturbations of for the case that the frequency map 0, is 
singular. 

2.1 Singularities of Lagrangian Maps 

In this section we briefly recall a few facts about the singularities of smooth maps |AGZV85] . 
ignoring the topology of the angle space. 

A map / : — > M" has a critical point at x if its Jacobian, Df{x), has less than maximal 
rank, i.e., if rank{Df) < min(m,?i). The image, f{x), of a critical point is a critical value. A 
map is said to be (C^) stable at x if every map that is sufficiently close to / (in the sense that 
the first k derivatives are close) is locally diffeomorphic to /. The equivalence class of maps 
that are locally diffeomorphic to / is the "germ" of /. If the dimension is low enough, the 
germ can be represented by a fixed polynomial map, more generally "moduli," which are either 
parameters or arbitrary functions, are needed to represent the germ. The equivalence class of 
maps represented by the germ at a critical point is called a "singularity." 

A Lagrangian map is defined by the projection of a Lagrangian manifold onto a Lagrangian 
plane. A prime example is geometrical optics where the Lagrangian manifold corresponds to a 
wave front together with its unit normals — the velocity vectors, and the projection is to physical 
space. Correspondingly, the set of critical values of a Lagrangian map is called a "caustic." 
A Lagrangian manifold can be represented by a single, generating function |Wei77j : if, as in 
our case, the Lagrangian manifold is a graph over the action space, {(p,J) = {DT{J),J), the 
generating function is T{J). The Lagrangian map is the projection of the manifold onto the 
action space, i.e. 9.:W^ defined hy J ^ DT{J). The map has a critical point at J if 

D^T{J) is singular. 

The standard theory of the singularities of Lagrangian maps has been formalized by Thom 
and generalized by Arnold p^GZV85 . When d = 1 there is only one stable singularity, the 
"fold." For d = 2 the "cusp" singularity is also stable. For d = 3, three new singularities are 
stable, the "swallowtail" and two forms of point singularities (pyramids and purses). 

2.2 Fold Singularity 

Suppose that d = 2 and that the generating function has the formal power series expansion 



so{^,n = {ip,i')+T{i') 



= DT{I) . 



(3) 




(4) 



i+j>i 



with the standard Euclidean scalar product (,). 
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The frequency map generated by T is denoted 

^1{J) = DT = {Ti,T2) = u* + 0{J) , 

where we denote the partial derivatives with subscripts i. Thus uj* denotes the image of the 
origin J = 0. The critical set is defined by the equation 

r(Ji, J2) = det{D^T) = T11T22 - Tf^ = . (5) 

The set of critical points is a smooth codimension-one submanifold in J-space, i.e. a curve, 
wherever the derivative Dt = (ti,T2) 7^ on the critical set. In this case, coordinates can be 
chosen so that the critical set goes through origin, and is tangent to the Ji axis, that is, so that 
^12(0) = T22(0) = 0. Under the nondegeneracy assumption that the Hessian D^T has rank one, 

Tii(O) = 2t2,o / O . 

In this case, the coordinates can be scaled so that ^2,0 = \^ ^-nd the generating function (jlj 
becomes 

T={u:\J) + \jl+ kjJl4- (6) 

i+j>2 

The singularity is locally a fold when the caustic is smooth — this occurs whenever the image 
of the tangent vector of the critical set is nonzero. A vector v is tangent to the critical set when 

Dt -v = Tivi + T2V2 = , (7) 
and its image, D^T ■ v, is nonzero if v is not in the kernel of D^T, 

ker(L»2r) = {w : Tnwi + T12W2 = , r = 0} . 

Therefore the function 

k(Ji, J2) = T11T2 - T12T1 

= 7ii(3Ti2Ti22 — ^11^222) + 712(^22^111 — 3ri2Tii2) , (8) 

is zero only when the tangent is in the kernel, and the caustic is locally smooth when 

«(J)|,(^)=o/0. 

This is the nondegeneracy condition for the fold. 

If we choose coordinates as in © , so that the singularity crosses the origin parallel to the Ji 
axis, then k(0) = — 24to,3. Thus there is a fold at the origin if the coefficient of is nonzero. 
After suitable coordinate transformations, all other cubic terms can be eliminated, and the 
coefficients can be scaled to give the generating function 

Tfoid = (^*, J) + ^iJf + J|) + 0(4). (9) 

Indeed, all higher order terms can be locally transformed away; hence, the generator © repre- 
sents the germ of the fold 'A (TZV85j . The fold singularity is denoted A2. 

The frequency map of © is Q{J) = ui* + (Ji, | Jf)) its critical set is the curve J*{s) = 
(s, 0), determined by det(-D^r) = 3J2 = 0, which is the horizontal axis. The caustic is the image 
of the critical set, = ^}{J*{s)), which a horizontal line. The action of the frequency map is 
to fold the J-plane into the upper-half 0-plane. 
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2.3 Cusp Singularity 

If K = at some point on the critical set, then it is tangent to ker(Z)^T). The tangency is of 
first order if, for any tangent vector ©, 

H{J) = Dk-v ^0 . (10) 

This is the nondegeneracy condition for the cusp singularity. 

In the coordinate system of ©, a cusp occurs at the origin when n{0) = — 24to,3 = 0. In 
this case, the cusp nondegeneracy condition ((TUl) reduces to 

/^(0) = 96ti,2(4to,4-t?,2)y^O. 

Consequently, the nondegeneracy conditions for a cusp are 

4to,4 tj ^ + 

Coordinates for Q can then be chosen so that all of the cubic and quartic terms apart from 
Ji JI and J 2 can be transformed away; collecting terms gives rise to the form 

T = (a;*, J) + (^Ji + ^ii,2J|) + \ (4*0,4 - t?,2) 4 + 0(5) 

After scaling, the cubic coefficient can be eliminated and the quartic one can be scaled to e = ±1, 
giving 

Tcusv = J) + \{Ji + Jif + eJl + 0(5) . (11) 

As for the fold case, the higher order terms are irrelevant. There are two cusp singularities, 
depending upon the sign e |A(tZV85] : they are denoted ^3. For (fTT|) . the critical set is the 
parabola J*{s) = (—(1 + 3e)s^, s), and the caustic is the semi-cubical parabola, or cusp: 




In the exterior of the cusp, the map is one-to-one, while in the interior it is three-to-one, see 
Fig. H and Fig. |2j The preimage of the caustic consists of both the critical set (shown in green 
in the figures) and the parabola Ji = — (1 + |e)J2 (shown in red). The tangent to the critical 
set is vertical at the origin, and it is mapped to zero at the cusp point. 

When e = 1 the critical set is mapped to the caustic by being, in effect, "rotated" about the 
negative Jli axis by 180°, while the other preimage of the critical set is simply squeezed, without 
rotation. 

For e = —1, however, the critical set is simply "squeezed" towards the positive axis, while 
the second preimage, a parabola opening to the left, is rotated by 180° and swept around to the 
positive caustic curve. 



3 Normal Form 

In this section we will show how to unfold the fold and cusp singularities — that is, obtain families 
of symplectic maps, 

{^',l') = F,4ip,I) , 
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Figure 1: The frequency map for the cusp ^4^ given by IIH with lu* — 0. The left figure shows the J-plane; 
the green parabola is the critical set, and the red is the second preimage of the caustic. The right figure shows 
the f2-plane; the caustic is the green cusped curve. The dashed fines show how a grid in J is mapped to O. 

on X that are perturbations of when has a fold or cusp singularity. To do this 
requires (at least) two parameters that we will call 5 and e. The first parameter will represent 
the unfolding of the singular set; it is the bifurcation parameter. We will assume that 5 = 
corresponds to the case that the singularity is at a distinct location, namely when the fold or 
cusp passes through a "resonance." The parameter e represents the strength of the perturbation 
that destroys integrability. When the dependence on the parameters is not essential we shall 
suppress it in the notation. 

The family F^^^ is generated by 

S{^,I') = {ipJ') + Ts{I')+eSi{^J',e,6) . (12) 

Here we will assume that the perturbation Si is a smooth function of its arguments, and that 
= DTs, has a fold or cusp singularity, thus, it has corank 1 on a critical set Ig{s) with curve 
parameter s. The caustic is denoted 

The perturbed map will in general also have a perturbed frequency map, determined by the 
independent portion of Si. However, for simplicity we may assume that all of the changes in 
the frequency map are encoded by the parameter S; specifically, we assume 

/ d^Si{ip,l',e,6)=0. 

Without loss of generality, the actions can be translated so that / = lies on the singular 
set, and in particular so that the origin corresponds to the point with the label s = 0; thus we 
can assume that /|(0) = 0, for any 6, see Fig. IHI The origin maps to the point = ^1(0) on 
the caustic. 

3.1 Unfolding from Resonance 

A frequency uj* is resonant if it is a rational vector: 

w* = (^,^), ni,n2GZ, iVGN, gcd(|ni|, jnaj, iV) = 1 . (13) 
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Figure 2: The frequency map for the cusp . The designation of the curves is the same as Fig. 

More generally, for a given a;*, define the sublattice 

£ = {m G : {m,LO*) G Z} . 

When dim£ = 0, then we say uj* is incommensurate, when dim£ = 1, lv* is commensurate, and 
when dimC = 2, oj* is resonant. The sublattice for a resonant example is shown in Fig. |1J 

We now argue that the parameters {s, S) can be selected so that uj* = loq is resonant. This 
is particularly important to construct a near identity approximation to the mapping, see ij3.41 
First consider the fold case. The parameter 5 unfolds the singularity if it causes the caustic 
to "move." For the fold this means that the map fi* : (s, 6) ^ uj from a neighborhood of zero 
in (s, 6) to the neighborhood of uj* in frequency plane is a local diffeomorphism. This ensures 
that the fold actually crosses uj* as 6 varies, recall Fig. 01 Since rational points are dense in 
the frequency plane, we can then define the labels s and 6, so that s = 6 = corresponds to a 
resonant 

UJ* = DToiO) = no{0) . (14) 

Note that we are still free to define the point s = on the singular sets for 5^0, and we will 
do so in the next section. 

We next select coordinates to align the angle variables with the "principle" resonance lines 
of UJ*. Let {fei,/c2} be a basis of the resonant module C, recall Fig. [l] Thus if we define the 
matrix 

A = {ki,k2), (15) 

then the sublattice is given by 

£ = {Ak:keZ^}. (16) 

There are a number of possible choices for these vectors; for example, they could correspond to 
the "minimal" basis — that with shortest length. They could also be chosen to correspond to 
the largest resonant Fourier coefficients in the perturbation, i.e., the directions of the "main" 
resonance lines that pass through uj* . In the next section we will argue that they might be 
selected to approximately align with the slope of the singular set. 

The matrix A defines a transformation to new angle variables (p = A^ip. A function U on the 
original torus T^, i.e. a function that obeys the periodicity condition U{if+m) = U[(p),\/m G 1?, 



Figure 3: Sketch of the (5-dependent frequency map fig for the case of a fold singularity. Here the origin in 
action space has been chosen so that / = is on the singular set and has constant slope. 

is transformed into a function 

u{^) = c/(A-V) 

on the new torus 

= m^/A^Z^ , (17) 

that is, a function with the periodicity condition U{(p + A^m) = U{(p),\/m G I?, see Fig. |21 In 
H',^A\ we will use this periodicity to average the angle dependence and eliminate nonresonant 
Fourier terms in Si. Generally the transformation A is not unimodular jCas97j : instead it has 
degree N = det A , meaning that copies of the new fundamental cell in (f are covering the 
original torus in ip. Equivalently, the unit cell in if has area jj. 

The angle transformation can be extended to a symplectic transformation, defining {ip, I) = 
(A^cp, A~^ I). As recalled in Appendix the generating function for the map in the new co- 
ordinates under this "point" symplectic transformation is simply S{<p,I) = S{A~^ip, AI), or 

S{p, T) = {p, T) + fs{i') + eSii^, /', e, 5) (18) 
where Ts{I) = Ts{AI) and Si{(p,I,e,6) = Si{A~^Lp,AI,e,6). Thus the new frequency is 

Co* = Dfo{0) = A^uj* . 

Suppose that x is the first order change in the caustic with 5, so that 

Co} = ns{0) = ^* + + 0{5^) . (19) 

The nondegeneracy condition ensures that x6 is not parallel to the fold — and freedom of choice of 
the point s = on the fold curves as 5 varies implies that only the component of x that is normal 
to the fold 17* (s) at s = is uniquely defined. Below we will select the unfolding parameter to 
be this component of x- For the moment, we retain the parameter 6 in the exposition for scaling 
purposes. 

For the cusp case, it is natural to shift the actions so that the origin of action space maps 
to the cusp point, at least for the special values s = 0, 5 = 0. Upon varying a single parameter 
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Figure 4: Sublattice L = Al? for tj = (|, |) with fci = (1, 2) and ki, = (—2, 3). These basis vectors were chosen 
so that fci • c<j* = gcd(|ni|, \ni\) — 1 and k2 ■ uj* = 0. (fci, te) is a minimal basis for the resonance module 

6, the cusp moves along a curve in frequency space that genericahy does not go through any 
rational point, thus the cusp frequency is not generally resonant. If however, 5 G and the 
map ^^1(0) is a local diffeomorphism from a neighborhood of (5 = to a neighborhood of uj* , 
then the origin s = 5 = can be selected so that the cusp frequency uj* is resonant. For this 
case both components of in (|19|) are relevant. In the final normal form, we will use these 
components as the unfolding parameters. 

3.2 Diagonalization of the Kinetic Energy 

~ (2) ~ ~ 

The second order terms in the kinetic energy correspond to a quadratic form (/) = (/, Msl), 
with "mass matrix" M^. By assumption has rank one, because the origin maps to the caustic. 
Let vs denote the unit vector along the kernel of M, 

M5V5 = , 

and Us the unit vector orthogonal to vg, i.e., the vector tangent to the singular set, recall Fig. 
Since Ms is symmetric, this implies Ms = XsusUg where is its nonzero eigenvalue. Thus the 
quadratic term in the kinetic energy can be written 

fPii) = Xs{usJf . 

The vector us, tangent to the singular curve at the origin, will rotate with 5. However, if the 
curvature of the singular set is nonzero at s = 0, the label s = along the singular curve can be 
selected so that us = uq. To see this, suppose that us{s) is the tangent vector to the singular 
set at s. Then, according to the implicit function theorem, there is a local solution s{6) to the 
equation us{s) x uq = precisely when ■§^us{s) is nonzero at the origin in {s,6). Relabeling 
the point s = on the singular sets for varying 5 makes u independent of 5. We call this the 
constant slope normalization. Note that this normalization also applies to the cusp case, where 
5 G M'^; however with this choice, the cusp is not necessarily the image of / = when 6^0. 
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Using the constant slope normalization, we can choose coordinates that are aligned with the 
kernel of Tj once and for all. Letting it* = (cos a, sin a) be the unit vector tangent to the 
singular set, we define 

(cos a — sina\ , , 

(20) 
sin a cos a j 

With this transformation the new kinetic energy becomes 

fs{i) = fsiRi) = {ui,Ri) + Xsil + ff\Ri) + . . . , 

so that the singular set is explicitly tangent to the Ji axis at the origin. 

To make the transformation (|2Ujl symplectic, it would be necessary also to transform the 
angle variables using (p = R^ip (recall Appendix^; however, this is not a map on the torus. 
To get around this problem, we transform the map with the nonsymplectic transformation 
((/3,/) = {ifjRI), leaving the angle variables unchanged. Under this linear transformation the 
generating function (fTH)) becomes S{(p,I') = S{(p,RI'), or explicitly 

SipJ') = {0,Ri')+fs{P) + eSi{ipJ',e,6) (21) 

Since the coordinates are not symplectic, the generating function cannot be simply differen- 
tiated to get the map. Instead, the map is now defined using (recall Appendix 1^ 

0' = R-'ArS, I = R"^-^S, 

or 

which yields the implicit mapping 

if' = ip + Cb* + 2A/ii?-*ei + i?"*V;,(f + ... + sSi) , 
i = P + eR-^V^Si . (22) 

Here the unit vector in the Ii direction is denoted by ei. This is a symplectic map written in 
non-symplectic coordinates. 
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If the slope tana were rational, then we could replace i? by a unimodular matrix whose 
second column is parallel to v to accomplish the effect that we desire. In this case we could also 
transform the angle variables without destroying periodicity. Note that if the null vector v is in 
the resonance module, then the same outcome could have been accomplished using the freedom 
in the selection of the matrix A, ()15|). thereby avoiding the transformation H2U() altogether.^ 

Even when the slope is not rational, we can select the module basis to nearly align with the 
null vector of M. In this case, the diagonalization can be accomplished with a rotation matrix R 
with a small. Thus the general case can be considered to be that of small a. Note however that 
the selection of the basis {ki,k2} of the resonance module with very large integer coefficients 
will mean that the Fourier terms with small integer coefficients in the new basis may correspond 
very large coefficients in the original basis. For "typical" functions, these coefficients will have 
small values and thus be less important dynamically. 



3.3 Scaling near Resonance 

To make explicit that we will study the neighborhood of the singularity at / = 0, we now choose 
a scaling so that I is small. The goal is to scale the action so that the terms in the kinetic energy 
that give rise to the singularity are of order e, the same as the perturbation. Naively, this would 
be attained by assuming that / = ©(y^), i.e., using the transformation ((^,1) = {'^jJ,^/eJ). 
Since this transformation is symplectic with multiplier (recall Appendix^, the new generating 
function is obtained by the transformation S{'ip, J) = e~2 S{'ip, \/eJ) or, from H21|) . 

J') = + Col RJ') + VI ^XsJ[^ + Siii;, 0,0, 6)\ +0{e) , 

showing that the quadratic term in the kinetic energy is now formally the same order as the 
perturbation. 

At this point the terms in T^^^ that give rise to the fold normal form Q are formally of higher 
order. Thus to "see" the fold, we modify the transformation by differentially scaling the actions 
in the directions tangent and normal to the singular set. Using the scaling transformation 

B = diag(l,z.-i) , (23) 

we redefine 

{^,i) = {^j,V^BJ) . (24) 

As before, the overall factor of y/e gives a transformation that is symplectic with multiplier, so 
that the new generating function is S{iIj, J) = -^S{xp, y/eBJ), or 

S{iIj,J') = {iP + ^l,RBJ') + ^\\5J[^ + Si{iP,0,0,5) \ +ei^~\s4^ + eOil,i^-\u-^) (25) 

When the coefficient k of J| is nonzero, the choice v = e^l^ makes the cubic kinetic term 0(-y/e) 
have the same order as the perturbation, and the neglected terms have higher order. Note that 
the nondegeneracy condition k 7^ is precisely that for the fold (jH)). 

To enforce the near resonance condition not only the action, but also the bifurcation param- 
eter xb, must be small. We can set b = \fe to make this explicit, leaving x to represent the 
unfolding parameter. It is also convenient at this point to introduce x using (|19j) . using the new 
representation 

^If the null vector of the mass matrix in the original coordinate system J is n G £, we would choose k2 — v, 
for then the transformed mass matrix M = A^MA has only one nonzero component. 
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so that the components of x correspond to the frequency shifts along and transverse to the fold, 
and the transverse component, X2 is scaled with J2. After this transformation we obtain 



5(V', J') = {i> + Cj\RBJ') + [(x, /) + Tfoid{J) + U{^)\ + 0{e^'^) (26) 

where C/(V') = 'S'i(V'i 0) 0, 0), and TfoidiJ) = AoJf + kqJI- Without loss of generality, we can 
scale the coefficients so that 

Tfoid{J) = l{J! + ^oJl) (27) 

by an overall choice of units for the action. The sign of the quadratic term can be changed by 
multiplying S and C by — 1. 

As was discussed in ^2.31 when k = 0, the singularity is generically a cusp. In this case, the 
same scaling transformation ()23() to select the leading order terms in the kinetic energy can be 
used, though with the choice v = e^l^ . There are now two additional terms that are of order 
\/e, one cubic and one quartic. The generating function again has the form (|26|1. though the 
kinetic energy becomes 

TcusviJ) = ^(-^1 + P^'^lf + , (28) 

and the first corrections are now Oic"!^). 

Since the additional transformation (|24|1 is also nonsymplectic, the map generated by S is 
obtained from the relations 

= , = cAs, (29) 



where 



^ 1 / cos a sma \ , , 

C^[^RB)-^= . . (30) 

\— z^sma i/cosa/ 



One clear advantage of the scaling is that the resulting map is now easily written explicitly or 
in the more compact leap-frog form 





V' + LD* + y^c* (x + vr(j')) 



(31) 



The complete set of transformations from the original variables (99, /) of (|12|) to the variables 
(-0, J) of (ED) is 



-0 = , 

J = ^CA-^I . (32) 
ye 



3.4 Iterating near Resonance 



The map F,.^^ (|31j) becomes a near-identity map after it is iterated N times, where N corresponds 
to the denominator of the resonant frequency w*, recall H13|). The point is that, all of the change 
in ()31() is ©(^/e), except for the term cD*. However, after N iterations, this term becomes Nuj* = 
A^Nlu* G J^1?\ therefore, this term is "integral" on the transformed torus and consequently 
disappears from the phase map after iteration. Since ^Je is assumed to be small, a near-identity 
map is obtained. 
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More precisely, using the evolution is 

i-i 

Jj = Jo-V^CY,^U{^l^t), 
Therefore ipj = V'o + 3^* + (^(-v/s); and the action after N steps can be simplified to 

N-l 

Jm = Jo- VIC ^U{iIjo + tCo*) + 0(e) . 
t=o 

This can be simplified by expanding the potential in a Fourier series. By assumption the 
perturbation is smooth, and thus so is the perturbing potential = 0, 0, 0). This 

function is periodic on the transformed torus T^, H17() : therefore it has the expansion 

Uii,) = c^e^-*^^"'™'^') . (33) 

Using the identity 



N 

2TTilp/N 



E 



\ip/N Z 
N else 



e 

1=0 

it is easy to see that only the resonant Fourier terms, those with m £ C, contribute to the 
Fourier series of Jn, so that 

Jn = Jo- VINCVU{iIjo) + 0(e) . 
Here U denotes the Fourier series of U with all the nonresonant terms removed: 

Uii;) ^ C(An)e^<"''^> . (34) 

Note that since all of the lower frequency terms have been removed from the Fourier series, U 
is now a function on the ordinary torus: 

+ k) = U{ip) , 'ike I? 

Thus, as far as the A^*^ iterate is concerned, we can take ijj mod 1. Consequently, the term 
Ntb* in the A^*^ iterate of the angle is equivalent to 0, and since Jj changes only at ©(-^e), we 
obtain 

^N = ^o + V^NC' ix + VT{Jn)) + 0(e) . 

We can conclude by noting that the near identity map for the N^^ iterate has the generating 
function 

Se,V^, J') = (V-, C-^J') + V^N ((x, J') + T{J') + i7(V)) + 0(e) . (35) 
The map is generated from H35|) using the nonsymplectic form ()3U() . 
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4 Hamiltonian Normal Form 



The main advantage of the near-resonance condition is that the A^-step map generated by H35|) 
is a near identity map; consequently, it is — to order e — the time N^/e map of a flow. This flow 
is generated by the Hamiltonian 

H{4^,J) = {x,J)+T{J) + U{i^). (36) 

Here T is either the fold, Tf^id (HZJ, or the cusp, T^ugp (|28|) ■ normal form. For the Hamiltonian, 
the nonsymplectic-form ()29() becomes a noncanonical Poisson bracket m. z = {ip,J), and the 
equations of motion are given by 



JVH , J = { (37) 




where is the Poisson matrix. The Hamiltonian flow preserves the energy, consequently the 
original map has an approximate conserved quantity when e is sufficiently small. The existence 
of a conserved quantity for a 4D symplectic map is nontrivial. 
The equations of motion for H36|) are 

iP = C\x + VT{J)) 
j = -CVU{ip) 

Using the form (|3U(1 for C, these equations can be written 

ipi = cos a(xi + dj-^T) - V sin a(x2 + ^jjT) 
^2 = sin a(xi + dj^T) + v cos a{x2 + dj^T) 
Ji = — cos ad^^ U — sin aS^j ^ 

J2 = y (sin ad^^ U — cos ad^^ U ) (38) 

It is tempting to introduce new angles V = but since the slope tana is generically not 

rational this transformation does not respect the periodicity of U, with the result that U would 
be quasiperiodic function of the new angles. 

The model ()38|) has several important parameters: the frequency mismatch vector x, the 
direction of the fold a, and the perturbation size 1^. In addition, the kinetic energy for the fold 
model depends upon the relative size of the fold term, kq and cusp model has the coefficients po 
and ^0- In addition there is of course the freedom of choice of the periodic potential U. Note 
that u is either e^/^, or e^/^ and is assumed small. If we take this to the extreme = then H38|l 
implies that J2 is a second constant of motion. Thus in some sense our system reduces to one 
degree of freedom. However, the angle dependence is quasiperiodic. We will discuss this case in 



5 Fold Model 

The system (jSHJ can be further simplified for the fold model where T = Tf^i^ is given by ()27() . 
In this case xi is irrelevant and can be set to zero without loss of generality by, for example, 
a shift along the Ji axis. The point is that xi corresponds to a frequency shift along the fold, 
but in the normal form the fold is a straight line. The parameter X2 is the essential detuning 
parameter. 
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If we assume that the original potential has strong Fourier components in the direction of 
the fundamental resonances, then the dominant terms in the potential will heU = —acos{^i) — 
bcos{ip2)-'^ In this case the Hamiltonian becomes explicitly 

Hfoidiip, J) = X2J2 + ^{Jf + KoJ|) - acos(V'i) - 6cos(V'2) • (39) 

Though this Hamiltonian appears to be separable into two individual one-degree-of-freedom 
systems, this is not the case when a 7^ because of the nonsymplectic form (jHZI)- 

It may be the case that additional terms in the resonance module will have amplitudes 
comparable to the ones included in (|39|1 . In this case, these terms should be included in U and 
the dynamics of the two degrees of freedom are coupled even when a = 0. 

The dynamics of the two-degree-of-freedom system (|39]) can be visualized using an appro- 
priate Poincare section. Though the standard sections are not complete, they can still give 
insights into the dynamics. For our purposes, it seems that the most useful section is the surface 
S = {ipi = 0, H = E}. This section can be projected onto the coordinates {ip2, J2) with a choice 
of sign for Ji, since the Hamiltonian (|39() depends quadratically on Ji. Two example phase 
portraits are shown in Fig. El 




Figure 6: Poincare section for the Hamiltonian ll^ with (i/'2,J2) G [— tt, tt] x [—3,2.5] and Ji > 0. The 
parameters are set to o = 6 = kq = 1, S = 2, a = O.OItt, and u = 0.05. For the left panel X2 = —2.5, and for the 
right X2 = -1.7 

In these sections there are two visible elliptic fixed points; they correspond to periodic orbits 
of the flow. For example in the right panel the orbit at {tp2,J2) = (0,-0.016) has motion 
in which Ji is oscillating in narrow region about 2.5 while ipi rapidly rotates. The orbit at 
0^2, J2) = (^)~1) has its first degree of freedom on an oscillatory trajectory about the origin, 
with Ji oscillating between ±1.4. For both of these orbits 1^2 and J2 remain nearly fixed. 

Our goal is to explain the main features of these phase portraits. 

5.1 Aligned Fold: a = 

The model (|89|) reduces to the one-dimensional case if the fold happens to be aligned with the 
resonant module, so that a = 0, for then the equations decouple. The dynamics in {tpi, Ji) 

^It is also possible that the original potential has no terms in the resonance module so that U — 0. In this 
case a higher order normal form must be used. This is the situation for some resonances (those with even period) 
in the standard nontwist mapping |Pet01irWCT04HWCT'M05| . 
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then correspond to a simple pendulum. In the decoupled {1^21 J2) dynamics, the parameter (if 
it is nonzero) can now be eliminated by rescaling time. The resulting canonical one-degree-of- 
freedom model, 

i?2(V'2, J2) = X2 J2 + y4 - 6C0S(V'2) , (40) 

is the same as that originally studied by Howard jHH95j and derived as a normal form by Simo 
|Sim98j . This system has equilibria at the points 

, „ ji±) , / 2x2 

For definiteness, suppose that ko,6 > 0. When X2 < there are four equilibria as shown in 
Fig. CJ-two are centers ((0, Jg"*^^) and (vr, ^)), and two are saddles ((0, ^) and (vr, Jg^-*)). 
These orbits correspond to the two families of resonant islands with rotation number lo* in the 
original coordinate system that are caused by the fold. These resonant islands are destroyed in 
saddle-center bifurcations at X2 = 0. The separatrices of the saddles undergo a reconnection as 
X2 passes through 

X* = -l{b'^o)'^ . (41) 

When X* < X2 < 0, there are "meandering" invariant circles, i.e., invariant circles that go above 
the upper island chain and below the lower one. 




-3-2-10123 -3-2-10123 



Figure 7: Phase space of the one-degree-of-freedom fold model, l|40|l for no = b — 1 and X2 =• The top line 
shows generic phase portraits for X2 ~ x*(l/2)^'^'', X*(3/2)'^'''', — X*(l/6)^'''^- The center phase portrait shows a 
meandering curve. The bottom line shows critical phase portraits for X2 = x*i 0. Contourlines shown are H2 — i/2 
for i = —5, —4, . . . , 5. 

The Poincare sections for the full Hamiltonian (|39() in Fig.lHJare similar to the phase portraits 
in Fig. [7| though since a 7^ in the former figure the two degrees of freedom are coupled and 
some orbits are chaotic. Since v <^ 1 the dynamics of J2 is much slower than that of Ji, and 
since a is also small, the coupling between the two is weak. It is interesting to note that the 
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chaotic orbits in Fig. IHl (in the regime of the meandering tori) correspond to values of {ipi, Ji) 
that are near the pendulum separatrix for this degree of freedom. Conversely near the center of 
the islands on the section the corresponding values of Ji are either large, so that the pendulum 
is in a rotating regime, or near zero so that the pendulum is deeply trapped. The island near 
'(/'2 = is created in saddle-center bifurcation that occurs X2 ~ —1.47 for the value a = O.Olvr 
shown. By contrast, the saddle-center bifurcation for the island near ip2 = still occurs near 
X2 = 0, as it does for a = 0. The reconnection bifurcation no longer occurs at a single parameter 
value as it does for the integrable model ()4U|) nevertheless, as can be seen in the left panel of 
Fig. a topological change in the the dynamics occurs just above X2 = —2.5 where the unstable 
manifolds of the saddle periodic points undergo a reconnection in their intersections. 

5.2 Nearly Aligned Fold: a small 

When V is nonzero but small, the dynamics of the the two degrees of freedom are coupled but 
J2 evolves slowly compared to Ji. In such a situation it is appropriate to apply an averaging 
method. This averaging is valid when the dynamics of the first degree of freedom are far from 
separatrices which would cause this motion to slow. 

A particularly simple situation occurs when a is small — indeed, as we argued in Hl-{.2l this is 
the general case. If formally a <C i^, then to lowest order H38|l decouples into a pendulum and 
the nontwist one-degree-of-freedom system (|4()|) . The phase portraits then will look similar to 
those of Fig. [3 just as we saw in Fig. IHl 

We can formally treat this case by treating a as a small parameter of the same order as v. 
Up to terms of first order in a, reduces to 

■01 = ^1 , 

3 

^2 = aJi + v{x2 + -J2) > 
Jl = —asvatpi — ahs\n.tjj2 , 

^2 = -1/6 sin 02 ■ (42) 

To lowest order we may neglect the small ip2 term in the equation for Ji, and the (0i, Ji) 
dynamics become a simple pendulum on the fast, 0(1), time scale. Thus this system has periodic 
orbits and has true angle-action coordinates. Denoting this true angle variable by then an 
average over can be performed whenever its rotation frequency is large compared to z^. The 
standard averaging theorem |LM88j implies that for time scales t ~ 0(i^~^) the slow averaged 
dynamics are given by 

3 

tp2 = a{Ji)^ + u{x2 + -^JI) , 

J2 = —h'bsm^2 ■ (43) 

where the average ()^ is performed over This is valid providing the first degree of freedom 
is not near the pendulum separatrix where the time scale separation breaks down. Barring this 
event, the eff'ect on the dynamics of (-02; J2) is a simple shift: instead of at X2 = 0, the resonant 
saddle-center bifurcation now occurs for 

X2 = --{Jl)'& 
V 

The critical value x* for reconnection will shift down from H41j) by the same amount. How- 
ever, it is important to realize that the value of J\ on the Poincare section is determined by 
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H{0, "02) J2) = E and so will vary with position, thus the bifurcation values corresponding to 
different orbits on the section will shift by different amounts. The simple prediction of a shift in 
the bifurcation values is verified by our numerical computations. In the example of Fig. IHl the 
elliptic island around (0, j'^^) corresponds to ipi in a rotating regime with Ji oscillating slightly 
about 2.4. Thus the predicted bifurcation point is X2 ~ —1.5, close to the observed value of 
— 1.47. By contrast the orbit around the elliptic island on the section at (vr, has Ji = 1.4 

which is in the oscillating regime for the first degree of freedom so that (Ji)^ = 0. Thus the 
predicted bifurcation value is zero, just as observed. 

The averaging theory fails when the pendulum motion in the first degree of freedom is near 
its separatrix. In this case the coupling between the degrees of freedom drives Ji to jump across 
the separatrix on the slow time scale. This causes slow chaos in the second degree of freedom 
and will be discussed more below. 




5.3 Quasiperiodic Pendulum: v = Q 

When a = 0(1) the two degrees of freedom of the flow H38|) are coupled in an essential way; 
however, since z/ <C 1 there is still a fast/slow spitting of the dynamics. In the extreme case that 
1/ = 0, J2 is constant and is a second invariant for the flow. Be that as it may, this is a singular 
limit: as we have seen earlier in this section, the dynamics with v nonzero can result in 0(1) 
changes in J2 on the time scale even when v is arbitrarily small. 

When = 0, (|38|) reduces to a Poisson system, z = JpV Hp{z), with the variables z = 
{ipi , ■02 J ) 5 the Hamiltonian 

Hp{'il)i,il)2, Ji) = ^Jf -acostpi - 6cosV'2 , 

and the Poisson matrix 

/ 

Jp= 

\— cos a - 

which is the upper left 3x3 block of the matrix J' in ()37() . This Poisson matrix has the Casimir 
invariant given by 

( ^« ) _ ( sma\ / Vi \ _ 

\i^vj y— sina cos a J \ip2 J 

consequently, the motion is along lines of slope tana in the angle space. When the slope is 
irrational, these lines are dense on the torus so the Casimir is only a local invariant. Along a 
line of constant the system has infinitely many center and saddle equilibria at the points 
where Ji = and 

5 tan a 

sm^/^i = sin02 • 

a 

These occur for a set of values of the energy Ep = Hp that is dense in one or more intervals in 
the range \Ep\ < \a\ + |6|. For example if a = 6 > 0, then there are centers for energy values in 
a dense subset of [—2a, 0] and saddles in a dense subset of [0, 2a]. When Ep is above this range, 
there are no saddles, and the motion is "untrapped." 

If we define a coordinate tpu along the line of the Casimir by ipu 

= cos a ■01 + sin a ip2, this 

system can be obtained from the one degree-of-freedom Hamiltonian 

Hqp{iPu, Ji) = -Ji — acos(cosa il^u — sin a ■0^) — 5cos(sinQ iIj^ + cosq -i/^t,) (44) 
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with canonical equations of motion in Ji). However, it is important to note that Hqp 
depends quasiperiodically on ipu unless tana is rational. We show typical contours of Hqp in 
Fig. El 



3 
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Vu 



Figure 8: Contours of the Hamiltonian 1441 for a = b — 1, tpv = and tana = 7, the golden mean. Shown are 
the energy levels of the seven critical points of Hqp for tpu in the interval [0, Stt]. 

5.4 Unaligned Fold: a > 0{u) 

The dynamics of (|39|) can be quite complex when a > 0{v). In this regime, application of 
averaging is problematic since the unperturbed system {u = 0) is the quasiperiodic pendulum 
H44() instead of the simple periodic pendulum. Multi- frequency averaging is valid only when the 
frequencies satisfy a Diophantine condition LM88. . Moreover, averaging fails for our system in 
the regime where the orbit passes near any of the infinitely many saddles of the quasiperiodic 
pendulum. Here we describe only some of the qualitative features of these dynamics. 

In Fig.EJwe show the effect of increasing a on the dynamics of ()38() . In these phase portraits, 
a ~ 3z/. The chaotic orbits in the central part of the section (near ^2 = 0) have dynamics that 
corresponds to multiple separatrix encounters of the quasiperiodic pendulum. There are also a 
number of island chains visible in this region that correspond to trapping in the potential wells 
of Hqp. For these islands the quasiperiodic phase il^u{t) is bounded but ranges over intervals 
larger than 27r. The large island near ip2 = corresponds to (V'i)'^i) trapped near the origin, 
and its dynamics are relatively simple because the v = {) motion of the first degree of freedom is 
then periodic. Similarly the tori for J2 < —1.7 correspond to values of Ji > 2 where averaging 
over the fast quasiperiodic pendulum motion is apparently still approximately valid. 

When E < Q the island near ^2 = tt is no longer on the section — it is energetically in- 
accessible, see Fig. E3 In this case the invariant, meandering torus on the section boundary 
corresponds to the first degree of freedom undergoing small oscillations about the origin. For 
negative energy, the orbits near the two fixed points (0, J^^) correspond to trapped motion in 
the first degree of freedom, and the chaotic portion of the section has moved to smaller values 
of J2. It is also clear from this figure that the saddle-center bifurcation of these two fixed points 
at '(/'2 = occurs at X2 = 0, just as in the model — a result consistent with averaging since 
{Ji)-d = for trapped orbits. 

As E decreases, the energy surface intersects the section 1/^2 = in two disconnected pieces. 
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Figure 9: Poincare section of the flow of 1391 for the same parameters as Fig. |S| except that a = O.OSvr and 
E — 0.5. For the left panel X2 = —1 and for the right X2 = —0.5 The vertical range corresponds to J2 £ [—2.2, 1.8]. 

One corresponds to the elliptic island around {^^J^^\ and the other to the nonresonant tori 
for J2 < J2 ^ ■ Since H39() is cubic in J2 , the energy can be arbitrarily small, but then the fold is 
not accessible. 




Figure 10: Poincare section of the flow of 1391 for the same parameters as Fig. |S| except that a = 0.057r, 
E = —0.5. For the left panel X2 ~ —1 and for the right X2 = The vertical range corresponds to J2 G [—2.2, 1.8]. 

When tan a is rational, the phase portrait exhibits resonances that correspond to this cou- 
pling. An example with tan a = ^ is shown in Fig. 111! Here a period-three island chain replaces 

the resonance near Jg"*"^ . The large degree of chaos in this section reflects orbits that repeatedly 
cross the separatrix of the {ipi, Ji) pendulum. As can be seen in the right panel of this figure, as 
X2 approaches 0, there are bifurcations creating a period-two island chain as well as reconnection 
bifurcations of the period-three islands, which finally are destroyed in a saddle-center bifurcation 
at X2 = 0. 
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Figure 11: Poincare section of the flow of H39^ for the same parameters as Fig.^lexcept that a = arctan(l/3), 
E — 1 and u = 0.01. For the left panel X2 ~ —1 and the vertical range corresponds to J2 G [—1.6, 2.0], while for 
the right panel X2 = —0.15 and J2 € [—1.4, 1.8] 



6 Cusp Model 



When the kinetic energy has a cusp singularity, T = Tcusp, (|28() and we choose the simplest 
potential U as before, the Hamiltonian (|36)) becomes 



The essential difference to the fold case is that for the cusp the equations do not separate even 
when a = 0; both derivatives of Tcusp depend on each of the Jj. Thus, even for the case a = 
and small i^, we must use averaging in order to reduce to one degree of freedom. 

On the fast timescale J2 is constant and the equations for {ipi , Ji ) are that of the simple 
pendulum with a shifted Ji. As for the case of the fold we can introduce the angle conjugate 
to the true pendulum action and average over it to obtain a slow system in (^2, J2) where Ji 
is simply replaced by its average (Ji)^. A necessary condition for Ji to be rotating is that the 
energy be sufficiently large. 

After scaling time to eliminate i^, the slow, averaged, one-degree-of-freedom Hamiltonian is 



Hcuspii^.J) = XiJi + X2J2 + T,{-h + PoJ^)'^ + ^J'oJ2 - acosV'i - 6cos^/;2 • 



(45) 



In this case, both components of x are important. 



6.1 Aligned Cusp: a = 

When a = the differential equations (|38j) become 



^1 = Xi + Ji + P0J2 > 

^2 = v{x2 + 2poJ2iJi + P0J2) + ^PoJl) 

Ji = — asin^i , 

J2 = —i'bsm'ip2 ■ 



(46) 




(47) 
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where f3 = 2/9o(Ji)^ and e = fiQ + Pq/2. Recall from that the two symplectically nonequiv- 
alent forms of the cusp are represented by the sign of e. Notice that negative e is only possible 
for /Uq < 0. By a scaling of J2, e can be reduced to ±1. The system has the discrete symmetry 

^^£,/3,X2(V'2, J2) = -iJ-e -/3,X2(V'2 + VT, -J2) 

SO that for the bifurcation analysis it is enough to consider, e.g., e = +1. The bifurcation diagram 
as a function of the two unfolding parameters X2 and /? for this case is shown in Fig. 1121 For 
e = — 1 the picture is reflected through the vertical axis so that the critical values for /3 all become 
non-negative. There are four bifurcation curves shown in Fig. 1121 The curve + 27x1 = 
with a cusp at the origin corresponds to a saddle-center bifurcation. The three remaining curves 
mark reconnection bifurcations.^ The lower line of phase portraits all correspond to X2 = 0; 
they are symmetric under J2 — > — ^2- The shown part of the bifurcation diagram with X2 ^ 
has four disjoint regions separated by curves of critical values. The lower half of the bifurcation 
diagram with X2 < is related to the upper half shown in Fig. El by the discrete symmetry 




Figure 12: Bifurcation diagram of the averaged one-degree-of- freedom Hamiltonian l)47|l with e = +1. The 
dots in the phase portraits mark isolated stable or bifurcating equilibrium points. There are four generic phase 
portraits inside the respective regions. The five codimension one phase portraits are touching the fat critical lines. 
The three codimension two points are at X2) = (0, 0), (-4^, 0), (-4^/2/3, (8/3)^''''8/9). 

Moving counter clockwise on a large circle around the origin starting in the 4th quadrant 
near /? = reproduces the twistless bifurcation associated to a fold: first the creation of two 
island chains in a saddle-center bifurcation and next the reconnection. The characteristic feature 
of the cusp appears close to the origin in parameter space where the lines corresponding to folds 
are close to each other. The three reconnections lines differ in which of the three different lines 
of equilibria is reconnecting. A qualitatively new phase portrait appears in the triangular region 
attached to the origin, which has a localized invariant twistless curve. The main feature of this 
new region is that there are twistless tori which are not rotational invariant circles between the 

^The reconnection curves are given by -4194304+8192/3'' -4/3* -1769472/3x2 + 1728/3^^X2-/3^1 - 186624/3^X2- 
81/3'^xt - 2187/3^X2 - 19683x1 = 0. 
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figure-eight separatrix and the rotational separatrix. Unhke the meandering curves (which are 
also present in the same phase portrait) these are, however, localized in ■02 • 

Each of these regimes is easily visible in numerical phase portraits if one chooses parameter 
values so that the motion of the first degree of freedom gives rise, through ( Ji)^, to an appropriate 
value for /3. In Figs. H^llTll we choose fiQ > and po < so that e = +1, and fix all of the 
other parameters of (|45() . The dynamics of the first degree of freedom is then determined by the 
value of the energy: when E is large, points on the Poincare section have large positive {Ji)^, 
which makes P large and negative. As E decreases, so does (Ji)^, moving the phase portrait 
horizontally in the bifurcation diagram Fig. 1121 Each of the four regimes of Fig. El are shown 
in Figs. I13I14I The degree of chaos in these panels also increases, to the final case where the 
dynamics of Ji undergoes separatrix crossings over a significant portion of the phase space. 




Figure 13: Poincare sections of the cusp Hamiltonian H45^ with parameters po — —0.5, /io — 0.875, a — b = 1, 
V = 0.05, a = 0, xi = and X2 = 0.5. For the left panel E = 35.75, and for the right E = 18.11. The vertical 
range is J2 G [—0.2, 1.0]. The energy determines the value of (Ji)^. 




Figure 14: Poincare sections of the cusp Hamiltonian 14511 for the same parameters as Fie. 1131 For the left panel 
E = 3.75, and for the right E = 0.75. The vertical range is ,h G [-1.7, 1.6]. 
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6.2 Nearly Aligned Cusp: Small a 

Treating a to be of the same order as the small parameter the first order approximation is 

-01 = Xi + >/i + P^A ) 

i>2 = i^{x2 + 2poJ2{Ji + P0J2) + ^pq4) + "(xi + ^1 + pqJI) , 

(48j 

Ji = — asin^i — ab sin ^2 , 
J2 = —i'bsmip2 ■ 

The Hamiltonian for these equations again is H45() , but the equations of motion are different from 
the a = case because of the non-diagonal non-standard symplectic form. Averaging works as 
before after neglecting the small ip2 term in the equation for Ji. After averaging the equations 
of motion, we see they are generated by the one-degree-of- freedom Hamiltonian in {ip2, ^2): 

^ = (X2 + ^(Xi + (^1)^)) J2 + Po{Ji)^JI + ^PoJI + {Po + pI/2)4 - &sinV'2 • 

Finally this Hamiltonian can be recast in the form H47() by completing the quartic, scaling, and 
collecting terms so that new effective parameters (3 and X2 are found. Denote the change of scale 
of J2 necessary in order to make the quartic coefficient equal to e = ±1 by s, = e/{po + Po/^)- 
Then the shift in J2 is A = —s^eapo/{l2i'). The shift in /5 = 2{Ji)^pos'^ is given by — 6eA^. 
Finally X2 is shifted by 

s (xi + W^^) + 2AsVo(Ji)^ - 8eA=^ . 

Since when ajv = 0(1), the shift A is small, the dominant term in these shifts is typically the 
one involving ( Ji)^. Thus in particular making a non-zero breaks the discrete symmetry present 
when X2 = 0. 

6.3 Quasiperiodic case: z/ = 

Setting 1/ = causes J2 to be a constant, and thus again we obtain a quasiperiodic system with 
one degree of freedom. The essential difference to the fold case (|44|) is that here there is a linear 
term in J\ in the Hamiltonian. Thus up to a shift of x\ + PJ2 ™ "^i figain the quasiperiodic 
pendulum is found. 

Because the critical values are dense in an interval the quasiperiodic pendulum also has an 
abundance of twistless curves. In any region of non-critical points in phase space bounded by 
the separatrices of two distinct hyperbolic points there is a twistless curve. 

7 Conclusion 

We have derived two normal forms for a nearly integrable, four-dimensional symplectic map with 
a singularity in its frequency ratio map, J 1— > $7(J). A singularity corresponds to the vanishing 
of the twist T = det(Dr2). For a four-dimensional mapping there are two stable singularities: the 
fold and cusp. In each case we focus upon a point 0. = uj* where the caustic crosses a resonance 
curve (m, f^(J)) G Z for m G 1? . We defined coordinates (Ji, J2) so that the Ji-axis is tangent 
and the J2-axis is normal to the singular set. Consequently, J2 represents the degenerate action 
and a singularity in the frequency map then corresponds to the vanishing of quadratic terms in 
the kinetic energy for J2. Using e to denote the size of the perturbation, the dominant balance 
along the singular set gives J\ ^ 0(e). If the singularity is a fold, then dominant balance 
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transverse to the singular set gives Jf ~ 0(e), while if it is a cusp then j| ~ 0{e). The 
differential scaling gives rise to the small parameter u = e^/^ or e^/^, respectively in our normal 
form. 

Because it is nearly resonant, the normal form has an approximate invariant — an energy — 
that restricts the motion to three-dimensional surfaces. The corresponding dynamics is given by 
the flow of a two-degree-of- freedom Hamiltonian system of the form H = {x, J) + T{J) + U{ip). 
As can be expected, this system is generally nonintegrable. One source of nonintegrability is 
standard: the potential U contains terms that couple the degrees of freedom. Indeed this occurs 
even in the case of nondegenerate frequency maps, where T(J) = ^{jf + |Pos93j . 

In our model, however, there are two additional sources of complexity. First, the cusp model 
has an intrinsic coupling between the longitudinal (Ji) and transverse (J2) actions because the 
kinetic energy T in (|28jl is not separable. The second source corresponds to the fact that the 
singular set is not generally aligned with the resonance module, and indeed will generically have 
an irrational slope. This gives rise to the misalignment angle a in and the corresponding 
nonsymplectic structure of our equations of motion H38() . Though it is common for issues of 
rationality to be important in Hamiltonian dynamics for rotation numbers, as far as we know, 
the rationality of a slope in action space has not been important in any other situation. 

Because of the disparate scaling in the action variables, the time scale for evolution of the 
degenerate action J2 and its conjugate angle are slower by the factor u than those for the 
nondegenerate ones. If we simply set u to zero, then J2 is constant. The dynamics for this case 
reduced to a system with one degree of freedom with a quasiperiodic potential. An equivalent 
description of this quasiperiodic pendulum uses a rank one Poisson-structure on M x T^. 

The simplest dynamics occurs for the fold when a = 0, when the fold curve is parallel to a 
vector in the resonance module, and when the coupling harmonics in U are negligible. For this 
case the multi-dimensional fold behaves as the well-known one-degree-of-freedom case derived 
by Simo |Sim98j . However, even the fold has complex behavior when a ^ 0. We have seen 
that the coupling of the nonsingular action Ji to the degenerate degree of freedom J2 leads to 
slow chaos when the first degree of freedom is near a separatrix. It also causes the shift in the 
reconnection bifurcation values when Ji is in a rotating regime. 

The cusp dynamics corresponds to the collision of two fold curves. Consequently, near a 
cusp there can be three sets of island chains with the same rotation number. As the cusp is 
approached there are several "reconnection" scenarios. The intrinsic coupling of the two degrees 
of freedom in this case make the study of this system more complex. Even when a = the 
two degree-of-freedom normal form flow is non-integrable. When v is small averaging allows 
for reduction to one degree of freedom. This simple system has been completely analyzed, and 
good agreement is found with the full equations for appropriate choice of parameters. A new 
type of twistless torus which is not a graph over the unperturbed torus appears near the cusp. 
Similar to the case of the fold the case of small a can again be reduced to the same one degree 
of freedom model, with shifted parameters, however. 

It would be interesting to study the non-integrable dynamics of the normal form flow with 
two degrees of freedom. It seems possible to still use averaging theory to partially reduce the 
dynamics, but the crossing of separatrices and the corresponding jumps in action would need to 
be taken into account. Another interesting case is the perturbation away from the quasiperiodic 
case 1/ = for arbitrary a using some kind of quasiperiodic averaging theory. 
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A Appendix: Nonsymplectic Coordinates 

A canonical generating function S{(p,I') implicitly defines a symplectic map by the equations 
This map preserves the symplectic form 

w = dl Adip , 

since w — w' = d{Idip — I'dip') = d{Idip + ip'dl') = d^S = 0. The symplectic form has the 
equivalent matrix representation 

/ -id\ 
id j ' 



w 



in (if, I) coordinates. 

Under a general, linear point transformation {ip, J) = {Bip, CI), we can define a new gener- 
ating function 

J') = 5(s-V,C"V') . 

The dynamically equivalent map in the new, nonsymplectic coordinates is 

i,- = BC±,S, J = CB>^S. 

Thus the transformation is symplectic if i? = C~*. 

Symplectic transformations with multiplier are a special case of the above where BC* is a 
multiple of the identity. 
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